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Functional  Relations  for  the  RandT  Operators  on  Plane-parallel  Media* 

by  Rudolph  W.  Prelsendorfer 
University  of  Calif omla^  La  Jolla,  California 
1*  Introduction  -  In  the  classical  theory  of  radiative  transfer,  one  may  as¬ 
sociate  with  a  separable  plane-parallel  medium,  l*et,  a  medium  In  which  the  phase 
function  jfl  m  4-9^ foe  is  It^ependent  of  depth,  two  fUnctlcwis  S  and  T, 
which  analytically  characterize  the  diffuse  reflectance  and  ttmnoilttance  of  the 
medium.  Si^llelt  knowledge  of  the  S  and  T  functions  Is  tantamount  to  a  solution 
of  the  nultlple  scattering  problem  associated  with  the  medium.  VIhereas  the  volume 
scattering  function  <r  and  the  volume  attenuation  function  oc  sunnarlee  the 
Local  optical  properties  of  the  medium,  the  S  and  T  functions  susmarlze  the  alcbal 
optical  properties  of  the  medium.  Furthermore,  as  the  looal  behavior  of  the  light 
field  Is  described  by  the  equation  of  transfer  In  terns  of  (T  and  ,  so  is  the 
global  behavior  of  the  U^t  field  described  by  the  principles  of  Invariance  In 
,'tems  of  S  and  T.  The  global  approach  was  initiated  by  Ambarsumlan,  but  It  remained 
'  .  for  Chandrasekhar  to  fonnulate  the  four  principles  of  invariance  and  deduce  from 

i  ..auf 

',them  the  requisite  Integral  equations  for  S  and  T  for  the  ease  of  finite,  separable 
T*  » plane-parallel  media.  The  inqportance  of  the  global  foneulatloo  has  been  lUustrated 
by  Its  role  In  many  successful  solutions  of  multiple  scattering  problems  which 
hitherto  were  Intractable  under  the  local  fonoulation. 

In  this  note  we  consider  the  case  of  finite,  ncn-senarable  plane-parallel 
media.  As  Is  shown  below,  the  study  of  these  media  requires  not  two  but  four 
functions:  a  reflectance-transmittance  pair  for  each  of  the  two  boundaries  of  the 
medium.  The  four  prinelpLes  of  invariance  for  a  non-separable  plane-parallel 
medium  are  formulated,  and  are  fouixl  to  be  Just  sufficient  to  yield  the  appropriate 
tntegro-dlfferentlal  equations  governing  the  four  funotlcms. 

^Contribution  from  the  Scrlpps  Institution  of  Oceanography,  New  Series  No. 
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Functional  Relations  for  the  R  and  T  Operators  on  Plane-Parallel  Media 

It  Is  of  some  interest  to  observe  that  the  following  derivations  can  be 
patterned  after  the  procedure  used  by  Chandrasekhar  for  the  separable  case,  namely 
the  procedure  which  makes  explicit  use  of  the  analytical  statements  of  the  principles 
of  invariance  in  conjunction  with  the  equation  of  transfer.  This  procedujre  promises 
to  be  adequate,  under  appropriate  extensions,  for  use  in  the  formulation  of  the 
functional  relations  governing  the  scattering  functions  in  general  one-parameter 
carrier  spaces.  In  this  way  a  useful  methodology  already  extant  in  the  classical 
theory  is  retained,  and  extended  to  more  general  settingsi  In  addition^  this 
procedure  appears  to  be  pertlnant  to  an  important  question  in  the  mathematical 
foundations  of  radiative  transfer  theory,  namely  the  question  of  the  apparent 
mathematical  ablogenesls  of  the  S  and  T  functions  and  the  principles  of  invariance. 
In  short,  are  the  principles  of  invariance  part  of  the  basic  laws  of  radiative 
transfer  theory  or  merely  consequences  of  them?  In  this  connection,  we  observe  that 
the  operator  form  of  the  equation  of  transfer  adopted  during  the  course  of  the 
present  derivations,  yields  a  pair  of  statements  which  clearly  represent  vdiat  might 
be  termed  the  local  forms  of  the  princijd.es  of  invariance.  These  in  turn  may 
possibly  supply  a  key,  by  means  of  the  general  theory  of  Green’s  functions  for 
linear  operators,  to  an  existance  proof  for  the  S  and  T  functions  and  perhaps  a 
derivation  of  the  statements  of  the  principles  of  invariance  themselves.  If 
this  or  any  coterminous  procedure  is  possible,  then  the  general  equation  of  transfer 
may  serve  as  the  sole  starting  point  of  an  unbroken  network  of  deductions  covering 
all  the  salient  mathematxcal  features  of  modem  radiative  transfer  theory. 

The  operators  used  in  the  present  note  are  left  in  undecomposed  form  in 
order  to  secure  mAxlrnmn  symmetry  and  compactness  of  the  formulations. 


3 


Functional  Relations  for  the  R  and  T  Operators  on  Plane-parallel  Media 

2,  Equation  of  Transfer  and  the  Local  Forms  of  the  Principles  of  Invariance  -  We 
begin  with  a  plane  parallel  medium  embedded  in  Boolidean 

Figure  x 


three  space,  namely  that  subset  contained  by  two  planes  at  ge^etrlc  depths  a 
and  b  from  some  datum  plane.  Several  intermediate  planes  at  depths  x,  y,  ■ 
vd.ll  also  be  considered;  a^x*y<  z^b.  The  volume  attenuation 
and  volume  scattering  functions,  •c  ,  and  CT  ,  along  with  the  radiance 
function  N  (specific  intensity)  will  depend  spatially  only  on  the  depth  parameters, 
and  in  an  arbitrary  way.  The  steady  state,  emission  free  equation  of  transfer  for 
radiance  may  be  vnritten 

where 

We  adopt  the  concepts  of  outward  and  inward  radiances:  K (-j 
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and  write  the  outward  and  inward  radiance  distributions  at  level  y  as  ,  W-fi)  j 

these  are  functions  vdiich  assign  to  the  pair  i  ^  ^  *  \  ,  o<  Zir  , 

the  appropriate  radiance  values  at  level  y  in  the  indicated  directions.  If  ^4. 
and  —  are  the  totalities  of  outward  and  inward  unit  vectors  at  each  level 


y  ,  the  equation  of  transfer  may  be  written  in  the  following  operator  form: 


, 

d'i 

where  for  each  y,  a  «.  y  3 


ttg)  N.d)  , 

^  -»•  nJ*  t*!) 

b,  and  arbitrary  ( ^  ,  56), 


(2) 


0<  ffc  i  2*  t 


77  ]  clfj'c/4'  ^ 


Thus  K3)  and  i  (<j)  evidentally  play  the  role  of  local  reflectance  and 

transmittance  operators,  CT-  and  0^4  are  appropriate  restrictions  of  CT  chosen 

so  that  their  ^-arguments  are  positive,  (Such  restrictions  are  possible  because 

of  the  implicitly  assumed  isotropy  of  the  medium).  The  pair  (2)  clearly  represents 

the  local  (or  integro-differental)  foms  of  the  principles  of  invariance, 

3.  Principles  of  Invariance  for  Non-separable  Plane  Parallel  Media  -  In  an  earlier 

note^  the  classical  principles  of  invariance  were  generalized  to  the  context  of 

general  one  parameter  carrier  spaces.  These  generalizations  were  based  directly 

on  considerations  of  the  classical  forms  as  given  by  Chandrasekhar,  and  resulted  in 

2 

an  explicit  analytical  embodiment  of  the  principle  of  Invariant  imbedding  for  the 


radiative  transfer  context.  The  two  main  statements  of  the  principles  of  invarieuace 
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were  shown  to  be  of  the  forms: 


I 

II 


(x^oc  $  2  <  b. 

(3) 


The  structural  similarity  of  the  sets  (2),  (3)  can  be  noted.  By  suitably  choosing 
the  spread  of  the  interval  [x,  z~]  and  the  location  of  y  within  it,  the  operator 
forms  of  four  principles  of  invariance  for  the  present  context  may  be  obtained 
from  (3) ! 


I  -  N+fb)Ta,y)  +  N/-<y)Rta,bK 

II  N.^J)  •»  N_(a)T^ajy)*  N4  ^3) *•)  . 

Ill  N-rb)  s  M.^A)  T(«.6)  +  Nf(b)  Rri>>a)  N-(3)T(y,b)  '*•  N4. ( b) R^b,y), 

IV  N+  (a)  *  N4  T  CbiO.)  4-  N-  ^0)  R  fa,b)  *  Nf  (b)  T  (g,a)  +  N-  ^o)  . 

For  every  pair  of  depths  (x,  z)  in  the  closed  interval  [a,  b^  we  have 

R(*.ai=  ^  Id^^'dr. 

Furthermore 

where 

*  X  [  3  d^'dd\ 

^  J  dy'd/', 

in  vdiich 

Tri*)  *  «=fp  f“-p' I  1  i  , 
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T(ir^  H)  is  the  transmittance  operator;  i)  ,  are  its 

reduced  and  diffuse  components*  R^x,2)  is  the  reflectance  operator*  and  since 
the  reflected  flux  is  by  definition  entirely  of  diffuse  character,  the  correspond¬ 
ing  reduced  component  R  *(  x,  i )  of  RCx/Z)  is  the  zero  operator ^so  that 

It  may  be  shovm  that,  in  analogy  to  the  S  auid  T  functions  of 
the  classical  theory,  the  present  R  and  T  functions  have  the  properties! 

RCX.Z,  =r  (r~l=C  -/4,4  , 


4*  SSi2ii2SSi  Relations  for  the  ^  and  T  operators  -  In  what  follows,  let  K-  (a) 
be  of  arbitrary  angular  structure  and  set  N4  (b)  w  O  (the  zero  function), 
furthermore,  all  functions  and  operators  will  be  assumed  continuously  differentiable 
with  respect  to  the  depth  parameter  y*  Now  from  I,  with  the  adopted  boundary 
conditions  in  force. 


dN-l'j) 


RO,h) 


N-  (':>) 


d  R(‘^,h) 


(4) 


Since  interest  is  at  present  centered  on  the  functional  relations  for  the  slab 
defined  by  the  planes  at  depths  a  and  b,  we  consider 


I  irm 


which,  by  (2),  has  the  form; 

-  +  N.  ^a)  1^(0)  3  =r  -  N- t>)  -t  -*■  j 

>diere  the  latter  relation  follows  form  use  of  the  left-hand  equality  of  III. 


Similarly 


im 


a-***- 


(\/.(»)  [  i(o  I  RCo.b)  . 
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Since  ,  and  is  arbitrary,  the  limit  operation  l,n^ 


3-**-  Ta 

applied  to  (4)  yields; 


3'»«- 


It 


=»  hro)  -f  it  (a)  -t-  R  (a, b)  if  a)  RCo,io) , 


Now  from  II, 

dN. 


'3 


rJ-ZQ)  drra.’jl  ^  ^  ,  c/ 


(5) 


By  (2)  and  the  left-hand  equality  of  IV, 


lina. 


irtn 


‘3'>t 


o(  Ni» 


-  N.foi  T f<3,b)  icb) 


N-  /-I  T<'6,t)  ^Cb)  . 


Applying  the  limit  operation  |,rv^-^  .  to  (5),  setting  dT(^ib)  c^Tfa,'^) 

-i  aT"  ^ 


and  recalling  that  is  arbitrary,  we  have 

in  3t (  0;b) 


at 


T(o,h)i(k)  +  T(o,t)  hfb)  Rfb.a). 


The  procedure  is  now  clear.  By  applying  this  procedure  in  turn  to  the  right-hand 
equalities  of  principles  III,  and  IV,  we  have  the  results; 

III,  9jU^)  ^  T(o,b)  y(b)  TCb,a)  . 

db 


IV>  -  aTfOlb) 
3a. 


R  (a,b)  yfa)T(a,b)  . 
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In  conclusion,  several  general  obseirvations  may  be  made,  (i)  For  non-separable 
spaces,  i.e»  spaces  in  which  the  phase  function  jfk  »  4Tr<r/o<.  is  dependent  on 
depth,  I'  •-  IV*  show  that  the  R  and  T  operators  possess  polarity;  R(a,b)  ^  R(h,a) 
^  T(bjA)  ^  system  I*  -  IV*  is  a  simultaneous  system 

of  the  simplest  kind.  It  can  be  solved,  in  principle,  by  successively  considering 
the  individual  relations  in  the  order:  I*,  IV',  III',  II'.  (iii)  To  obtain  the 
equations  governing  the  R  and  T  functions,  one  must  first  assume,  as  in  the  classical 
case,  a  Dirac-delta  structure  for  hi-ia)  :  fj.ca)  =»  M*  •  Then 

by  successively  applying  the  operators  on  each  side  of  the  statements  I'  —  IV* 
to  and  making  the  appropriate  reductions,  the  requisite  integiv-differential 

equations  for  the  R  and  T  functions  are  obtained. 
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